FRW and Bianchi type I cosmology of f-essence 
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Abstract F-cssence is a generalization of the usual 
Dirac model with the nonstandard kinetic term. In this 
paper, we introduce a new model of spinor cosmology 
containing both Ricci scalar and the non minimally cou- 
pled spinor fields in its action. We have investigated the 
cosmology with both isotropy and anisotropy, where the 
equations of motion of FRW and Bianchi type-I space- 
times have been derived and solved numerically. Finally 
the quantization of these models through Wheeler-De 
Witt (WD) wave function has been discussed. 

Keywords Spinor Cosmology, Spinor Fields, FRW 
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1 Introduction 

Astrophysical data of supernovae of type la indicate 
that we l ive in an accelerated expansion era of th e 
Universe (|Perlmutter S.I Il999t iRiess A.G. et all Il998l ) . 
There are two major theoretical explanations for this 
phenomena: the first category is the fluid model in 
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which we keep the Einstein gravity as a dominant the- 
ory and introduce a fluid in the right hand side of the 
Einstein field equation and the n investigate the cos 



mological e volution of the mode l (jCopeland E. J. et al 



2006h and (| Amendola L.II2OO0L paper II). Another ap- 
proach which has been investigated in recent years is 
the geometrical one in which we search for generalized 
models of the gravity whic h can deduce the acceler ated 



expansion of the Universe (jElizalde E. et al 
per I). Some examples of later model is 
itv dNoiiri S. et al. 112003.2007.20111 ICognola G et al 



20081: lAzadi M. et all 



200 



2004 pa- 
f(R) grav- 



( Elizalde E. et al. 



2004, 



pape r III), |jamil M. et aL^Oll , paper I) and ( Momeni D. et al"^ 



2009L paper I) in which the only dynamical sector of 
action is a function of the Ricci scalar. Here the non- 
linear terms of the curvature can be regarded as an 
alternative for the accelerated expansion of the Uni- 
verse. Another class of models are the f(R,G) mod- 
els, where both curvature and Gauss-Bonnet terms as 



the dyna mical quantities (lElizalde E. et al. I I2004L pa- 



20 10:) and ( Felice A. et al. 



per I I) , (jBamba K. et al 

2010h . Another interesting approach is the f(T) grav- 



ity where T is the torsion. Although it has no cur- 
vature, the space-time manifold has time evolution, 



only 


(Moller C. et a 


. 1961. 1978: Havashi K. et al. 


1967.1973.1977.1979 


Wu P. et al.l20] 


0.201ll Pellegrini C. et al. k 


1963 


: lYerzhanov K.K. et al.l 20ld 


Chen S. H. et al.1 


2011 


: iBeneochea G. R. et al. 12011: 


Li B. et al. 2011: 


Wei H. et al. 


2011). Since its equations of motions are 



lower order, working with it is easier, also the matter 
comes from a non-minimally action and curvature. This 
model can be explained as the trace of the specified 
energy-momentum tensor, the common R i cci sc alar, 
and the Lagrangian f(R , T) lHarko T. etldl (|201lh and 



( Mo meni D. et al. 2009. Daper II). If we limit to the 
scalar field theories, there are many o ptions to have 



an accelerated universe as quintessence (jZlatev I. et al 
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1999; ISteinhardt P. J. et al.lll999) and ( Amendola L.I Variation of Lagrangian (5) with respect to a 
200C , paper I) quintom models d Feng B. et al.. I 2006 ; the equation of motion o f the scale factor 
Zhao W. kooalNoiiri fc Odintsov Il2006allbl ). K-essense jArmendariz-Picon C. et al.| 
1999.2000.200l|l and other combined models. There 
is no limitation for the usage of the fermions in 
the matter sector of the theory as a matter source. 
There are some works on spinor cosmology in the 



2ad - 



afK = 0. 



yields 



(6) 



literature (lArmendariz-Picon C. et al 
last paper 



Ribas M. O. et al 



( Vakili B. et al 
20051 ^ 



1999.2000.2001 



The variation of Lagrangian (5) with respect to ■0, ip 
is the corresponding Eulcr-Lagrangian equations for the 
fcrmionic fields 



20051 paper I I ) and 



Chaichcrdsakul K.I 12007 



Boehmer C. G.I 120081) . In this work, following an ear- 



lier work on the gener alization of the spinor cosmology 
( Jamil M. et al.ll201ll paper II), we have introduced a 
new model of spinor cosmology, in which the action 
contains both R and the non minimally coupled spinor 
fields. We also have investigated both isotropic and 
anisotropic cosmological models in this framework. We 
derive the basic equations of motion and solve them for 
FRW and Bianchi type-I models numerically. Finally 
both classical and quantum models in this framework 
have been discussed and some analytical solutions for 
the quantum cosmology have been studied. 



2 FRW metric in f-essence 

In this section we would like to present the derivation 
of the equations of motion for FRW metric in the f- 
essence. 

Let us consider the following action of f-essence 



5' 



(1) 



where R is the scalar curvature, Y is the kinetic term 
for the fermionic field ?p and K is some function of its 
arguments. In the case of the FRW metric 



ds 2 



-dt 2 + a 2 {dx 2 + dy 2 + dz 2 ), 



R and Y have the form 



R = 6 



(2) 



(3) 



1.5-AVyV 
a 



0. 



+0.5AV7V 

AV^7° 
+O.5AV07 
Another equivalence form is 



1.5-AVV7 
a 



iKi 



0. 



(7) 



(8) 



3#i^ + 



0. 



4- K$ = 



K, 4 -, = 0. 



Also the zero-energy condition is given by 
Lad + L^ip + L-ip - L = 0, 
which yields the constraint 
YK Y - K = 0. 



-3a _2 d 2 



(9) 



(10) 



(11) 



(12) 



Collecting all equations and rewriting using the Hubble 
parameter H = (In a) t, we obtain a system of equations 
of f-essence (for the FRW metric case): 



3ir 

2H 



-p = 0, 
3H 2 +p = 0, 



2,HK 4 , + K,J: + K^ 



+K 



SHKj + Kj.tjj 



p + 3H(p+p) = 0. 



0. 



0. 



(13) 
(14) 

(15) 

(16) 
(17) 



Y = O.5i(V>7°0- V^V), 



(4) 



respectively. Substituting these expressions into (1) 
and integrating over the spatial dimensions, we are 
led to an effective Lagrangian in the mini-superspace 



L = -2(3aa 2 - a 3 K). 



(5) 



Here 



YK\ 



AT, p = K, 



(18) 



are the energy density and the pressure of f-essence. It 
is clear that these expressions for the energy density and 
the pressure represent the components of the energy- 
momentum tensor of f-cssence as: 



Too = YK Y - K, T u = T 22 = T 33 = —K. (19) 
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3 FRW model with K(Y, i/>, ij>) = Y - V(il>t/>) 

We introduce a useful model which is more applicable 
and more suitable for exact solutions. This model is 
described by 



K(Y,i/},ip) = Y- V(ip^). 



(20) 



To obtain the field equations, we substitute this form 
in Eqs. (13- 17). Another simple method is re-deriving 
these equations using the action directly, therefore in 
each of these equivalence methods we have the following 
equations of motion for Dirac fields: 

3H(lijrf) + - \ih° ~ = 0, (21) 



3ff(-±i7»ty) = V$. 
FRW equations in this case are 

2- + (-) 2 + kW-fyfy) 

a a l 

-W) - o, 

3H 2 = V(^). 



(22) 



(23) 
(24) 



The general potential is V(ipi>) = 2tpip. For this special 
case, we have the next set of EOMs: 

3H(-itp^°) + ^ij°i> - -i^7° — 2-0 = 0, (25) 



2H + -i(^7°V> - 07V) = 0. 



(26) 
(27) 



Now we take the Dirac 2-spinor as ip = (tpi, 1P2) 7°, the 
equation for spinor reads as 

i> + ^Htp + 2i 7 V = 0. (28) 
Thus we must solve the next system of ODEs: 

^ = |^ a ±2 l a = {l,2} = { +> -} > (29) 

which posses the following solution 

i? = (M0)a(t)^e 2t ,M0)a(t)^e- 2t ). (30) 

Using this form of the 2-spinor we can obtain the scale 
factor from the following equation 



M 2 e 4t }, 



(31) 



here y = bg(o(t)), {a,/3} = {V>i(0), ^a(0)}. There is 
no simple analytic solution for y(t). But if we take 
a = P = ^5, then we can solve it numerically. FIG.l 
shows the time evolution of y(t) for some initial values. 



4 Bianchi type I cosmology of f-essence 

The action of f-essence reads as 

5 = J d 4 x<S=g[R + 2K(Y,i>,i;)], (32) 

where K is some function of its arguments, ip = 
(ipl, ip2, ip3, i>A) T is a fermionic function and i> = ■0 + 7 ° 
is its adjoint function. Here 



Y = 0.5#r"I>^ - {D^)T^} 



(33) 



is the canonical kinetic term for the fermionic field and 



Dp is covariant derivative 



Dai) = dai) + ftui), Dai) = dui) - V^l 



(34) 



Here are spin connections, are the Dirac matri- 
ces associated with the space-time metric satisfying the 
Clifford algebra 



(35) 



The T M are related to the flat Dirac matrices, 7°, 
through the tetrads e° as 



= e ^ a r 



e u7a- 



(36) 



At the same time, the spin connections fi^, satisfy the 
relation 



f!, = 0.25 fa A(3,^ + ^<)/7 o . 



(37) 



The tetrads can be easily obtained from their definition, 
that is 



(38) 



Let us now consider the Bianchi type I universe filled 
with f-essence. These models for the special simple 
spino rs have been discussed previously (|Saha B. et al 



20041) . The metric is given by 

ds 2 = -N 2 (t)dt 2 +a 2 (t)dx 2 +b 2 (t)dy 2 +c 2 {t)dz 2 , (39) 

where a(t), b(t), c(t) are scale factors in the x, y, z direc- 
tions respectively and N(i) is the lapse function. The 
corresponding scalar curvature takes the form 



R = 



2 d b c 
iV^a + b + c 



ac 
ac 



ab 
ab 

be aN bN 



cN. 



be aN bN cN 



(40) 



where a dot represents differentiation with respect to t. 
For the metric (39) the tetrads take the form 



diag(N, a, 6, c), 

= diag(l/N, 1/a, 1/6, l/c). 



(41) 
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These formulas yield 



^ 01 r> _ ^ ..0..2 



° 0= °' ° 1= ~2]V 77 ' " 2 = ~27V 77 
2N ' 



(42) 



where 7 and 7* are the Dirac matrices in Minkowski 
spacetime and we have adopted the following represen- 
tation 



[ i r l a 1 



(43) 



Substituting (40) and (42) in (43) and integrating over 
the spatial dimensions, we are led to an effective La- 
grangian in the mini-superspace {N, a, b, c, ip, ip} 



L = -2[—{abc + abb + abc) - NabcK(Y, ip, ip)] , (44) 
where 



where energy density and pressure take the form 

p = N 2 (YK Y - K), p = N 2 K. (53) 

The vacuum solutions of the above system i.e. ip = 
is the generalized Kasncr solution in which space is 
homogeneous and has Euclidean me tric depending on 
time according to the Kasner metric ( Kasner. Edward] 
19211 ). This solution possess ing a Belinsky-Khalatnikov -1 



Belinskv V.A. et al.lll969l 



Lifshitz (BKL) singularit y 
Lifshitz E. M. et al.lll960h . which is a model of dynamic 



evolution of the Universe near the initial singularity 
t = and described by an anisotropic homogeneous 
and chaotic solution to the Einstein's field equations of 
gravitation. The Mixmaster universe exhibits similar 
properties as the Kasner solution. 

Some properties of g-essence were studied in 



Kulnazarov I. et"aL^ 



201lURazina O.V. et al. II2011I ). Model (1) admits two 
important reductions: k-essence and f-essence. 



L f = 2NabcK(Y,^,i>), Y = — (^V " ^V)- 



(45) 



The preliminary set-up for writing the equations of mo- 
tion is now complete. 

Variation of Lagrangian (44) with respect to TV, a, b, c, 
and ip yields the equations of motion of the gravita- 
tional and the fcrmions fields as: 



5 Bianchi type I cosmologies with 
K(Y,ip,$) — Y — V{ipip) 

In this section we examine the Bianchi type I cosmology 
for K (Y, ip, ip) = Y — V{ipip). In this case we have 

Vl ab ac be 



ab ac be 
ab ac be 
-N 2 (YK Y - K) = 0, 

b c be 
b c be 

N ,b c N , 

a c ac 
a c ac 

N a c 
b a ba 
b a ba 

N ,b a, n 
-N^b + a^ NK -^ 
K Y ip + 0.5[(]n(abc)) t K Y + K Y ]ip 
+N 1 °K^ = 0, 

K Y P\) + 0.5[(\n(abc)) t K Y + K Y ]Pp 
-NK^ a = 0, 

P + :iH(p + P ) = 0, 



(46) 
(47) 
(48) 

(49) 
(50) 

(51) 
(52) 



ab ac be 
-N 2 V(i>ip) = 0, 

b c be 
b c be 

-§(!;+ w^=°< 

N b c 

a c ac N ,a c . 
- + - + — --(- + -) 

a c ac N a c 

+N 2 K = 0, 

b a ba N b a 
b + a + ba~~ N^b + a> 
+N 2 K = 0, 

tp + 0.5(ln(a&c))tV> 
-N 7 % = 0, 

^ + 0.5(ln(a6c)) t ^ 

+NV^° = 0, 

p + 3H(p + p) =0. 



(54) 

(55) 

(56) 

(57) 
(58) 

(59) 
(60) 
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We take Vfyt/j) = 2t/jijj. Thus we obtain 

ab ac be 
ab ac be 
-2N 2 $i> = 0, 

b c be N ,b c i 
o c be -/V o c 



1 



+ N ^ ~ ^ ^ " 2?w = °' 

ci c ac N ,a c N 

- + - + T7 (" + -) 

a c oc ly a c 

b a ba N ,b a „ 

7 + - + ]T" + 
o a oa iv o a 

+N 2 (^$-y°i> - h°i>) - 2^) = o, 

i) + 0.5(ln(a6c)) t V - 2A/W = 0, 
^ + 0.5(ln(a6c))t^ + 2N^° = 0, 
p + 3H(p + p) =0. 



(62) 



(63) 



(64) 

(65) 
(66) 
(67) 



We solved these equations numerically for a set of ini- 
tial conditions imposed on the set of the functions 
{a(t),b(t),c(t), \^(t)\ 2 }. We set N(t) = 1 without loss 
of generality, since the metric is a projectable metric, 
i.e. we can define a new time coordinate t' = J N(t)dt. 
The gauge N(t) = 1 chosen in classical cosmological 
models, and called the cosmic time gauge. Another 
gauge fixing leads to N% =0, here iVj is the shift vec- 
tor. 

The numerical solutions is shown in the FIG. 2. As 
we observe, the functions {a(t), b(t),c(t)} are monoton- 
ically increasing functions of time, but the density func- 
tion IV^i)! 2 is a decreasing function of t. 



6 Quantization of f-essence 

For quantization of the model as descr ibed in (32) we 



adop t the method proposed by Misner (j Misner C. W. 
19691) . The first step is writing the general Hamiltonian 
suitable for describing the quantum evolution of the 
system. For a typical model (32), with an unknown 
form of the function K, it is not possible to write such 
Hamiltonian. But if we restrict to the case K(Y, ip, %p) — 
Y — V(iptp), the problem at hand becomes tractable. 
Th is is a special case of t he form discussed previously 
in (jVakili B. et al. 1120051 paper I). First we introduce 



a set of the metric functions 

x= 27! log( ^ 



( 61 ) p = log(\/ 'cVab), 



1, ,Vab 

° = o lo 8'( ) 

2 c 



(68) 

(69) 
(70) 



It is easy to show that the Hamiltonian for model 

K(Y,^,i>) = Y-V(tpip) is 



H 



,-3p 

~V2 



( P 2 p -pi-pi) + e i P[V(^)-T}=0, (71) 



where T is the Lagrange multiplier of the system. 
The set of the corresponding conjugate momentums 
of the new set of configurational coordinates {p, er, x} 
is {pp,Pa,Px} which satisfy the commutation brack- 
ets. With a specified form of the interaction Vtyip), 
we can obtain the classical solutions described in 



( Vakili B. et al. 1120051 paper I). Now from (71), we 
get the wave functio n directly from the WD equation 
( DeWitt B. S. Ill9671 ). and with the usual replacements 
Pi — > ~ «-g§-, the Wheeler-De Witt equation is 



r 1 
12 



{-d 2 + dl + d 2 ) + e^(-e- 3 " - T) *(p, a, X ) = 0. (72) 



We write the wave function as ^(p, a, x) = e t ^ k " a+kxX ^r(p)^ 
where 

- d 2 T(p) + [-hi - k 2 x + 12e 6 "(-e- 3 " - T)]T(p) = 0. (73) 
The general solution for (73) is 
-1/3 iV5 



r(p) 



-3/2 p 



d M 



-C«W 




1/3 ik, 4/3iv / 3v / Te 3p (74) 



, 1/3 ifc, 4/3iv / 3v / Te 3p 



Here M and W are WhittakerM and W hittakerW 
functions respectively while k - 
total wave function is 



kl + k\. Thus the 



V(p,a, X ) = J2 e t(k " a+kxX) e- 3/2p 



(75) 



a k M( l J^H } i/ 3ifc) 4 /3 iv / 3VTe 3p 



+ b k W 



-1/3 i\/3 



, 1/3 ifc, 4/3iv / 3v / Te 3p 



For normalization we set bk = and truncate the series 
for convergence. Thus the solution is written in final 



6 



form as 
*(P.0-,X) = 



(76) 



1/3 ik, 4/3iV3Vre 3p 



The Fourier- Whittaker coefBcients can be obtained 
from the initial wave function ^(O, a, x)- For Gaussian 
wave packet we can obtain the following result for the 
Fourier amplitude as follow: 



2 71 



(77) 



In Eq. ([77| we use the approximation T<1. 



7 Conclusion 

In conclusion, we derived the equations of motion of 
f-essence for FRW and Bianchi type I metrics. It is 
shown that if the Lagrangian of fermionic fields K has 
the usual Dirac form than the corresponding results 
coincide with the standard Einstehi-Dirac theory. We 
have investigated both classical and quantum aspects 
of this model. 
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Appendix 

As a double check, one can obtain the above field equa- 
tions from the Einstein and Dirac equations given by: 



= T, 



TOD^ + Kq = 0, 
D^T^ + K^ = 0, 
p + 3H(p + p) = 0. 



(78) 
(79) 
(80) 
(81) 



For a homogeneous fermionic field ip(t), equations (A2) 
and (A3) are equivalent to (Al) and (A4) respectively. 
On the other hand, the non-vanishing components of 



the Einstein tensor for the metric (39) are 
Goo 



G*22 

G33 



ab 


ac 


be 




ab 


ac 


h ^ 




a 2 


-b 


c 


be 


N 2 


.6 H 


h 

c 


be 


b 2 


a 


c 


ac 


N 2 


.a 


h - + 

c 


ac 


c 2 


-b 


a 


ba 


N 2 


lb" 


h 

a 


ba 



N ,b 
~N b ^ 
N /a 



N fb 
N\b 



(82) 
(83) 
(84) 
(85) 



The components of the energy-momentum tensor for 
the fermionic field as the matter source can be obtained 
from the standard definition as: 



dL 



f 



Q g y,v 



(86) 



yielding 



T Qa = -2N 2 (YK Y ~ K), T 11 = -2a 2 K, T 22 = 
T 33 = -2c 2 K, Ty = T 0i = 0. 

Substituting these results into Einstein equations (Al), 
yields the same equations as (46) -(52). In the case 
of the FRW metric (3), the equations corresponding to 
the action (1) can be obtained as: 



-2b 2 K, 
(87) 





3H 2 - p 


= 0, 


(88) 




2H + 3H 2 +p 


= 0, 


(89) 


Kyi) + 0.5(3HK Y - 




= 0, 


(90) 


Kyi) + 0.5(3HKy - 


VKy^ + lK^ 


= 0, 


(91) 




p + 3H{p + p) 


= 0, 


(92) 



where the kinetic terms, the energy density and the 
pressure take the forms 



Y = 0.5i(ipj°ip-ipj°ip), 



p = K Y Y - K, p = K. 



(93) 



(94) 
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'^£7 1 Numerical solution for y(t) = log(a(i)) . 
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Fig. 2 Numerical solution for {a(t), b(t), c(t), \i>{t)\ 2 } 



